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OS . Abstract 

^ . In bistable systems, the stability of front structures often influences the dynamics of extended pat- 

CX) ' terns. We show how the combined effect of an instability to curvature modulations and proximity to a 

pitchfork front bifurcation leads to spontaneous nucleation of spiral waves along the front. This effect 
is demonstrated by direct simulations of a FitzHugh-Nagumo (FHN) model and by simulations of order 
parameter equations for the front velocity and curvature. Spontaneous spiral-wave nucleation often re- 
sults in a state of spatio-temporal disorder involving repeated events of spiral wave nucleation, domain 
spliting and spiral wave annihilation. 
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Figure 1: The evolution of an unstable front solution connecting the up state (black) to the down state 
(white) in a bistable FHN reaction-diffusion system. The frames, from left to right, represent the pattern 
solution at successive moments in time. Perturbations on the flat traveling front grow and nucleate spiral 
pairs along the front line (see the Appendix for a different view of the middle frame) . The resulting disordered 
state is characterized by domain splitting and spiral-wave nucleation and annihilation. 

1 Introduction 

Spatio-temporal disorder is a ubiquitous phenomenon in extended nonequilibrium systems yet the mecha- 
nisms responsible for it are only partially explored. A generic mechanism for the onset of disorder in periodic 
patterns consists of phase instabilities followed by the formation of phase singularities which appear as dis- 
location defects or vortices [1]. This mechanism has been observed in numerical simulations of the complex 
Ginzburg-Landau equation [2, 3], and in experiments, e.g. electroconvection in liquid crystals [4]. A con- 
siderable effort has been devoted to elucidating the nature of the transition from the regular phase regime 
to the regime where vortices or defects spontaneously appear [5, 6, 7, 8, 9]. Other mechanisms involving 
instabilities of periodic patterns, leading to spiral breakup, have been reported recently [10, 11, 12, 13]. 

In this article we present a mechanism for the onset of spatio-temporal disorder associated with front 
structures. This type of disorder is illustrated by a numerical solution of a bistable FitzHugh Nagumo (FHN) 
type reaction-diffusion system (Fig. 1). An almost flat front, connecting the two stable uniform stationary 
states, begins traveling through the system. The front represents an "up" state (black) invading a "down" 
state (white). Initial nonuniform perturbations of the front position grow into wiggles which nucleate pairs 
of spiral waves (see the Appendix for a different view of the middle frame in Fig. 1.) The solution then 
evolves into a disordered state with repeated events of domain splitting and spiral-wave nucleation and 
annihilation. Similar phenomena have been observed in experiments on the ferrocyanide-iodate-sulfite (FIS) 
reaction [14, 15]. 

2 The NIB front bifurcation and spontaneous front reversals 

The key to understanding the disorder associated with front structures is a pitchfork bifurcation in the 
velocity of a propagating front. The front bifurcation is represented by the equation 

C 3 -(A c -A)C = 0, (1) 

where C is the velocity of a flat front and A is a control parameter. The corresponding bifurcation diagram 
is shown in Fig. 2. A stationary front (C = 0) becomes unstable below a critical parameter value, A = A c . 
At that point two new stable front solutions appear, representing an up state invading a down state (C > 0) 
and a down state invading an up state (C < 0). Hereafter, we refer to these front solutions as to "UD front" 
and "DU front" , respectively. 

This type of front bifurcation has been derived for periodically forced oscillatory systems [16] and for 
bistable FHN type models [17, 18, 19]. Experimental observations of front bifurcations, or supporting 
evidence for their existence, have been reported in Refs. [20, 21] for liquid crystals, and in Refs. [22, 23] 
for chemical reactions. The stationary and counter-propagating fronts are sometimes referred to as Ising 
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Figure 2: The pitchfork front bifurcation. At A = A c the stationary (C = 0) front solution becomes unstable 
to a pair of counterpropagating travelling fronts. 

and Bloch fronts, respectively, and the bifurcation itself as a nonequilibrium Ising-Bloch (NIB) bifurcation 
[16, 18]. 

Physical realizations of front bifurcations usually involve perturbations that unfold the pitchfork form of 
(1) into [24] 

C 3 - (A c - X)C + ui + v 2 C 2 = . (2) 

An asymmetry between the up and down states is an example of such a perturbation. In the following we 
confine ourselves to the case ^2 = 0. A plot of the surface (2) in the space spanned by C, A and u\ is shown 
in Fig. 3. The significance of small variations of v\ in the vicinity of the pitchfork bifurcation point, A = A c , 
v\ = 0, is now evident: perturbations may induce transitions between the upper and lower sheets and reverse 
the direction of front propagation. Notice that farther from the bifurcation point the variations of v\ must 
be larger in order to induce front reversal. 

The dynamics of a single flat front may also be changed by slow dynamical processes such as an approach 
to a boundary, an interaction with another front, or development of curvature. Fig. 4 shows a typical graph 
of the normal velocity of a front, C„, versus its curvature, n, for fixed A near the front bifurcation. The figure 
was obtained using a singular perturbation analysis of an FHN model, assuming slow curvature dynamics 
with respect to the time scale of front reversal [25] . The hysteretic shape, similar to the graph of C versus v\ 
for fixed A (see Fig. 3), suggests that front reversals can be induced by small curvature perturbations. Since 
curvature is an intrinsic dynamical variable, reversals of propagation direction occur spontaneously as the 
curvature of a front changes. Spontaneous front reversals in catalytic reactions on platinum surfaces have 
indeed been observed for parameters in the vicinity a front bifurcation [22]. Experimental observations of 
front reversals induced by boundaries have been reported in [23]. 

Imagine now a flat UD front that is unstable to transverse perturbations (i.e. an instability to curvature 
perturbations). As the front propagates, alternating segments along the front acquire negative and positive 
curvatures. For parameters that place the system near the the front bifurcation, where Fig. 4 applies, 
segments with negative curvature will eventually reverse propagation direction. Such local front reversals 
involve the creation of transition zones between the counterpropagating UD and DU fronts. These zones 
form the cores of rotating spiral waves. 
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Figure 3: Top: The surface C 3 — (A c — A)C + v\ = in the space spanned by C, A, v\ (the cusp catastrophe). 
First triad: A section at v\ = showing the pitchfork bifurcation (center), and sections at v\ < (left) 
and v\ > (right) showing unfoldings of the pitchfork. Second triad: Sections of the surface at constant A 
showing the hysteresis point (center) with single valued (left) and multivalued (right) relations away from 
the hysteresis point. 
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Figure 4: Normal front velocity C n vs curvature k in the vicnity of the front bifurcation. The development 
of small negative curvature may induce a transition from a UD front (C„ > 0) to a DU front (C n < 0). 
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Figure 5: A closer look at front interactions in the numerical solution of the FHN model from Fig. 1. The 
repulsive interaction between approaching fronts causes them to reverse direction. The reversal is followed 
by domain splitting. 

The scenario sketched above provides a heuristic explanation for the spontaneous nucleation of spiral 
waves in Fig. 1 which precedes the onset of spatio-temporal disorder. Similar arguments hold for other 
intrinsic perturbations such as front interactions. Indeed Fig. 1 includes many events where local front 
reversals are induced by the interactions between approaching domains. An example of such an event is 
shown in Fig. 5. Similar processes have been observed in experiments (Figure 6) on the FIS reaction [14, 15]. 



3 The dynamics of front reversals 

Algebraic relations like that displayed in Fig. 4 are useful for predicting the onset of spiral-wave nucleation: 
nucleation events become feasible when a single valued relation becomes multivalued (or hysteretic). Such 
relations do not, however, contain information about the nucleation process itself. To study the nucleation 
process a time-dependent approach for the front dynamics is needed. Near the bifurcation, the asymptotic 
front dynamics are governed by both a translational degree of freedom and an order parameter associated 
with the bifurcation: the front velocity C. In Refs. [26, 27] we derived asymptotic dynamic equations for a 
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Figure 6: Patterns in the Ferrocyanide-Iodate-Sulfite reaction show interactions leading to front reversals 
followed by domain splittings. These images are from experiments performed in the Center for Nonlinear 
Dynamics at the University of Texas at Austin. 



single unperturbed one-dimensional front using FHN type models. The equations for the front position, X, 
and velocity are 

X = C, (3) 
C = (a c - a)C - /3C 3 . 

These uncoupled equations describe the convergence to constant speed motion Ising (zero speed) and Bloch 
fronts. Order parameter equations for front propagation in one space dimension near a NIB bifurcation have 
also been derived in Ref. [28] to study the effects of external fixed heterogeneities. 

The effect of curvature is to couple the two equations in a way that allows for front reversal. The equations 
for the dynamics of a single two-dimensional front with smooth curvature, k, are [29] 

and 

^ = K-a)C-/3C 3 + 7 K + 7o + ^-^^ KC n ds' , (5) 

where C„, the normal front velocity, is given by 

C u = C-Dk, (6) 

and s is the arclength coordinate along the front. In deriving these equations an asymmetry between the 
up and down states has been introduced. Equation (4), for the curvature of the front, follows from purely 
geometric considerations [30, 31]. Equation (5), for the speed of the front, is valid near the front bifurcation 
and the boundary of instability to transverse perturbations [25, 29]. The integral term in both equations 
represents "advection" of changes in C and k from the stretching of the arclength over time. Note that away 
from the front bifurcation where the time scale associated with front reversal, (a c — a) , is short, C is no 
longer a slow variable and can be eliminated adiabatically. For a circular front, equation (5) then reduces to 
an algebraic relation between the normal front velocity and its curvature such as the one in Fig. 4. 

We have computed numerical solutions of equations (4) and (5) starting with an almost flat UD front 
as an initial condition. Fig. 7 pertains to parameter values in the Bloch regime where the UD and DU 
fronts coexist and are both stable to transverse perturbations. The initial front converges to a flat UD front 
propagating at constant speed. Crossing the transverse instability boundary causes perturbations on the 
front to grow. Near the NIB bifurcation, the growing curvature triggers the "nucleation" of a front segment 
with opposite velocity as shown in Fig. 8. The front structures in the C-s plane, separating segments with 
positive and negative velocities, pertain to spiral waves in the physical two-dimensional plane. 
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Figure 7: A solution to equations (4) and (5) when the two-dimensional front is stable to transverse pertur- 
bations. The frames in the top row show the time evolution of C(s) (thin line) and k(s) (thick line). The 
solution converges to a constant speed flat (/c = 0) traveling front. The frames in the bottom row display 
the corresponding dynamics of the front in the physical two-dimensional plane. 
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Figure 8: A solution to equations (4) and (5) when the two-dimensional front is unstable to transverse 
perturbations. The frames in the top row show the evolution of C(s) (thin line) and k(s) (thick line). The 
frames in the bottom row display the corresponding dynamics of the front line in the physical two-dimensional 
plane. A small curvature perturbation grows and the negative curvature triggers the formation of a region 
where the front speed becomes negative. The boundary points of this region form the cores of new rotating 
spiral waves. 
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4 Conclusion 



We have presented a new mechanism for spontaneous spiral- wave nucleation in bistable media that leads to 
spatio-temporal disorder. Unlike most other mechanisms which involve destabilization of periodic patterns 
(see however [32, 33]), this mechanism involves destabilization of fronts and may induce spatio-temporal 
disorder from a single front state. The dynamical equations, (4) and (5), for the front speed and curvature, 
describe the asymptotic behavior of fronts near the front bifurcation and the transverse instabilities. These 
equations capture the process of spiral-wave nucleation and can be used to analyze the transition from the 
stable curvature dynamics shown in Fig. 7 to dynamics involving nucleation events shown in Fig. 8. It 
would be interesting to find if an intermediate parameter range exists where the curvature fluctuates but no 
nucleation events occur. 

5 Appendix 

The model used to demonstrate the nucleation of spiral-vortices in Fig. 1 is a doubly-diffusive version of the 
FitzHugh-Nagumo Equations 

u t = e^{u - u 3 - v) + 5~ x u xx , 
v t = u- a x v -a + v xx . 

The parameters a and a\ are chosen such that the equations have two stable uniform solutions (bistable) 
and e and S control the type and stability of front solutions between those two stable states. 

Fig. 9 shows a closeup of the middle frame in Fig. 1 The transverse instability of the original front 
solution has already caused the formation of spiral- vortex pairs along the front. The spiral- vortices are 
identified with the crossing points of the zero contour lines of the u and v fields. At these points the normal 
front velocity is zero. On either side of the spiral-vortex the front propagates in opposite directions. 
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